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摘 ”要 : 针对 飞机 结 冰 的 特殊 凝固 过 程 ， 基 于 
法 ， 并 基于 凝固 实验 开展 了 所 建 模型 及 预测 方法 的 验证 。 研 究 表明 ， 相 对 于 传统 Enthalpy-porosity 方 法 ， 所 发 展 的 方 
法 可 有 效 表征 非 平衡 凝固 过 程 ， 因 而 对 冷水 滴 凝 固 特性 的 表征 具有 较 好 的 改进 。 基 于 所 建 方法 ， 开 展 了 水 滴 凝固 特性 


的 影响 因素 分 析 ， 获 得 了 不 同 过 冷 条 件 下 水 滴 凝 固 过 程 的 温度 分 布 及 相 界 面 变化 特征 。 相 关 砂 


Enthalpy-porosity 模 型 ， 发 展 了 改进 的 结 冰 特 怕 


(中 国 空气 动力 研究 与 发 展 中 心 空 气动 力学 国家 重点 实验 室 ， 绵 阳 621000) 


: 163443 


EFH 


I 模型 及 数值 计算 方 


的 改进 ， 以 及 结 冰 特性 的 精细 化 预测 提供 重要 参考 。 
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究 可 为 结 冰 热 力学 模型 


Abstract: Based on the Enthalpy-porosity model, an Improved icing prediction model and numerical method 


are developed for the special freezing process of aircraft icing. The model and prediction method are 


validated by the related experiments. Compared with the traditional Enthalpy-porosity method, the results 


show that the new method can characterize the non-equilibrium freezing process of supercooled water 


droplet more effectively. Using the developed method, the influencing factors of the freezing characteristics 
of supercooled water droplet are analyzed, the temperature distribution and phase interface characteristics 


are obtained under different cooling conditions..The results can provide an important reference for the 


improvement of icing thermodynamic model and finer prediction of icing accretion. 
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结 冰 物理 过 程 的 复杂 性 ， 因 而 目前 大 多 关于 过 冷水 


滴 结 冰 的 研究 把 重点 放 在 凝固 第 二 阶段 ， 即 相 界 面 


飞机 结 冰 是 过 冷水 滴 撞击 飞机 机 体 表 面 3 


发 生 ”的 推进 过 程 上 %1。 鉴 于 此 ， 大 多 研究 者 基于 


凝固 的 特殊 液 / 固 相 变 现象 ， 也 是 影响 飞行 安全 的 重 固 的 相关 理论 与 方法 开展 结 冰 特 性 的 预测 研究 ， 


的 阶段 ， 即 校 晶 形成 阶段 ， 也 有 研究 者 称 其 为 部 分 司 过 程 的 影响 特征 。 


"d 


凝固 阶段 四 。 第 二 阶段 为 完全 凝固 阶段 ， 即 液 固 相 


界面 由 液 相 相 向 回 相 推进 直至 凝固 完成 的 过 和 


程 ， 使 凝固 的 第 一 阶段 显著 快 于 第 二 阶段 G。 


由 于 前 非 平衡 凝固 特性 预测 手段 


来 越 多 的 关注 953。 但 由 于 凝 


要 隐患 之 一 03]。 在 飞机 结 冰 过 程 中 ， 过 冷水 滴 的 凝 不 少 研究 者 8B, 9 基于 Enthalpy-porosity 法 研究 凝固 ; 
固 通常 分 为 两 个 典型 阶段 由 。 第 一 阶段 由 形 核 开 程 的 液 固 相 变 行 为 ， 但 该 方法 主要 针对 平衡 凝固 
台 ， 是 水 滴 从 热力 学 非 平 衡 态 过 渡 到 热力 学 平衡 态 


程 ， 无 法 表征 过 冷水 滴 凝 固 的 第 一 阶段 即 非 平衡 ; 


ct 
zn 


究 还 较为 缺乏 ， 
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Se nm 
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近年 来 ， 随 着 结 冰 预 测 精度 要 求 的 提高 ， 结 冰 
301, 过 程 的 过 冷水 滴 的 非 平衡 凝 
在 飞机 结 冰 过 程 中 ， 由 于 撞击 引起 的 异 相形 核 过 


司 效应 引起 了 研究 者 越 
固 过 程 的 复杂 性 


目 
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测 精度 也 有 待 于 提高 ， 因 而 过 冷水 滴 凝 固 规律 特征 
的 相关 研究 也 较为 薄弱 。 本 文正 是 针对 飞机 结 冰 过 
程 过 冷水 滴 凝 固 的 非 平 衡 效 应 ， 发 展 了 过 冷水 滴 凝 
固 特 性 预测 方法 ， 并 自行 搭建 了 实验 系统 开展 了 所 
建 方法 的 实验 验证 。 相 关 研 究 可 为 结 冰 热力 学 模型 
的 改进 ， 以 及 结 冰 特性 的 精细 化 预测 提供 参考 。 
1 非 平衡 凝固 模型 

实验 研究 表明 ， 过 冷水 滴 的 凝固 过 程 由 两 个 典 
型 阶段 构成 。 在 第 一 阶段 即 枝 唱 形成 阶段 ， 过 冷水 
滴 温度 由 过 冷 态 上 升 至 凝固 点 ， 并 伴随 潜 热 的 部 分 
释放 ， 是 非 平衡 凝固 阶段 ， 第 二 阶段 为 由 热 扩 散 驱 
动 的 界面 推进 阶段 ， 也 是 在 等 温 下 进行 的 平衡 凝固 
过 程 。 尽 管 研究 表明 凝固 第 一 阶段 相对 于 第 二 阶段 
的 时 间 相 对 较 短 ， 但 由 于 结 冰 条 件 的 不 同 ， 使 第 一 
阶段 结束 时 形成 了 第 二 阶段 的 不 同 初始 条 件 ， 从 而 
影响 了 凝固 的 后 续 特 征 。 因 此 ， 过 冷水 滴 整 个 凝固 
过 程 的 有 效 预测 应 综合 考虑 第 一 阶段 的 非 平衡 凝固 
效应 及 第 二 阶段 的 平衡 凝固 效应 。 针 对 过 冷水 滴 的 
凝固 两 个 典型 阶段 的 特点 ， 如 何在 相 界 面 推 进 过 程 
的 预测 中 考虑 第 一 阶段 的 非 平 衡 效 应 ， 是 有 效 预测 
结 冰 全 过 程 凝 固 特性 需 解决 的 重要 问题 。 


图 1 ”过 冷水 滴 凝 固 的 两 个 典型 阶段 


Fig. 1 Two typical freezing stages of supercooled water droplet 

考虑 到 凝固 第 一 阶段 完成 并 形成 混合 态 的 过 程 
中 ， 过 冷水 滴 已 经 有 了 潜 热 的 部 分 释放 ， 因 此 ， 凝 
固 第 二 阶段 预测 中 应 进行 考虑 凝固 第 一 阶段 特性 参 
数 的 修正 。 将 无 量 纲 过 冷 度 & 表示 为 : 


e=cAT/L 


(D 
Ap, c 为 固 相 比热容 ， 工 为 固 / 液 相 变 潜 热 ，AT 
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的 计算 表达 式 为 : 

Afer £o d seno 

(2) 

其 中 ，7 为 液 相 凝 固 温度 ， 丰 ， 为 过 冷 态 温 
度 


凝固 第 一 阶段 结束 时 混合 态 中 未 发 生 相 变 的 液 
相 分 数 矿 可 表示 为 6 


pata; (3) 


其 中 ，c 为 固 相 比热容 。 


因此 ， 混 合 态 的 液 / 固 相 变 潜 热 Le 可 表示 为 : 


Ls -L-.f 
(4) 
混合 态 的 有 关 物 性 参数 PM mi 可 表示 为 : 


PM „x = PM,f, +PM l- f) 


(53 
其 中 ，PM ,和 PM .分别 表示 对 应 的 液 相 和 固 相 物 
性 参数 。 


当 获 得 凝固 第 一 阶段 凝固 过 程 的 物性 参数 后 ， 
考虑 过 冷 阶段 非 平衡 凝固 效应 的 相 变 问 题 即 可 转化 
为 平衡 凝固 问题 。 为 借鉴 平衡 凝固 条 件 下 液 / 固 相 变 
传 热 的 预测 方法 ， 作 以 下 几 点 假设 : 

(1) 忽略 水 滴 凝 固 过 程 水 滴 的 变形 ; 

(2) ”由 于 撞击 形成 的 异 相形 核 作 用 ， 凝 固 第 一 
阶段 时 间 显 著 小 于 第 二 阶段 ; 

(3) ”凝固 过 程 的 液 相 区 实质 为 已 在 凝固 第 一 阶 
段 发 生 了 部 分 相 变 的 混合 态 ， 因 此 凝固 第 二 阶段 计 
算 过 程 液 相 区 物性 参数 由 混合 态 参数 代 蔡 。 

基于 上 述 假设 并 借鉴 Enthalpy-porosity 方法 ， 


则 凝固 过 程 的 流动 及 传 热 控制 方程 可 描述 为 : 
Qh 

Ox, 

(6) 


WORD 批 量 转 PDF 工 具 -未 注册 
注册 码 购买 联系 QQ:3049816538 


Ou, Ou; Op O | Ou, 
p +u, - +u 
Ot Ox, Ox Ox, | Ox; 


-pg,B(T -T,)- S, 
CD 


以 直径 为 4mm、 高 度 为 2.5mm 的 过 冷水 滴 为 
例 ， 分 析 其 凝固 过 程 的 温度 分 布 及 相 区 分 布 特性 。 


其 中 ， 水 滴 底 面 为 冷却 面 。 
示 。 


» 


计算 中 底部 为 等 温 边 
流 边 界 。 计 算 物 性 参数 如 表 1 所 


表 1 ， 冰 和 水 的 材料 物性 参数 品 


Table 1 


Oh oT 0 oT 
P Eum C Mn P NC 
ôt ' àx, ôx, Ox, 


(8) 
Kb, us u Mixo x, 分 别 表示 不 同方 向 的 速度 


> 
LÉ 


坐标 (对 于 三 维 问题 ，i, 7 =1,2,3 )。p 和 了 分 


别 对 应 于 压力 和 温度 变量 ， 


T PISA. hA. o. u. BASATA 
对 应 于 密度 、 粘 性 系数 、 热 膨胀 系数 和 热 导 率 等 物 
ESZ g, 为 不 同方 向 的 重力 加 速度 ， 为 凝固 第 


-— 


二 阶段 两 相 区 的 源 项 ， 可 描述 为 
S, 2 -C(1- f, Ju, 
(9) 


结 冰 过 程 凝固 第 二 阶段 液 相 率 的 比例 ， 可 表示 为 : 
0 7T<T 
T-T 
= : I SIXTI 
fh Der s L 
1 T>T 


(10) 


凝固 时 的 温度 。 
采用 有 限 体 积 法 离散 控制 方程 ， 非 稳 态 项 采用 
时 间 一 阶 精度 格式 ， 采 用 强 隐 式 迭代 法 进行 求解 。 
基于 固定 网 格 技 术 ， 不 直 


免 直接 追踪 两 相 界 面 带 来 的 困难 。 


基于 上 述 方法 ， 


了 方法 的 改进 并 开展 了 3 
2 非 平衡 凝固 模型 


型 工 况 的 算 例 分 析 。 


T, 为 满足 Boussinesq 近 


其 中 ，C 为 表征 凝固 过 程 两 相 区 特征 的 参数 ， 广 为 


Physical properties of water and iceD] 

HE pikgm? ci kg-K! A/Wm KA Lk kg! — 
水 . 1000 
UK — 920 


4187 0.6 
2031 2.3 


334624 


图 2 显示 了 过 冷水 滴 凝 固 过 程 的 温度 分 布 特 


MES 图 3 显示 了 过 冷水 滴 凝 固 过 程 的 相 区 分 布 特 
性 。 在 凝固 初始 阶段 ， 相 界面 与 等 温 线 与 底部 冷 表 
面 近似 平行 。 随 着 凝固 过 程 的 推进 ， 相 界面 与 等 温 
线 逐 步 向 水 滴 顶 部 弯曲 。 


0s 4s 8s 
12s 16s 30s 
S0s 70s 90s 


图 2 考虑 过 冷 效应 的 水 滴 凝 固 过程 温 度 分布 特 性 
Fig.2 Temperature distributions of droplet freezing 


considering supercooling effect 


其 中 ， 了 和 工分 别 为 液 / 固 相 变 过 程 完 全 熔化 和 完全 


i 接 追 踪 两 相 界面 ， 通 过 温 
度 与 丛 值 的 关系 ， 反 算 不 同时 刻 的 相 界 面 位 置 ， 直 


对 项 目 组 前 期 自行 编制 的 相 变 传 热 计算 代码 0 进行 
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图 3 考虑 过 冷 效应 的 水 滴 凝 固 过 程 相 区 分 布 特性 


Fig.3 Phase interface distributions of droplet freezing 


considering supercooling effect 


图 4 显示 了 过 冷 度 为 -10C 条 件 下 考虑 与 不 考虑 


非 平衡 凝固 效应 的 凝固 速率 比较 ， 图 5 显示 了 水 滴 


不 同 过 冷 度 对 凝固 相 界 下 


变化 速率 的 影响 特征 。 可 


以 看 出 ， 在 考虑 过 冷 期 间 非 平衡 凝固 效应 影响 条 件 


下 ， 凝 固 速 率 较 传统 方法 有 所 提高 。 这 是 由 于 所 建 
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模型 考虑 了 过 冷 阶段 形 核 和 校 晶 生长 过 程 中 水 滴 潜 
热 的 部 分 释放 现象 ， 因 而 在 界面 推进 阶段 ， 液 / 固 相 
变 潜 热 的 释放 量 相 对 减少 ， 因 而 界面 推进 速度 较 传 
统 方法 有 所 提高 。 从 图 5 可 以 看 出 ， 在 相同 条 件 

下 ， 水 滴 过 冷 度 越 大 ， 相 变 的 驱动 力也 相对 越 大 ， 

因而 凝固 速率 也 相对 越 高 。 


1.2 


E P 
= y pet 
P 一 
€ » A 
8 2- 
tz ¥ w^ 
$ 0.4 E 
E y vj 
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图 4 过 冷 效应 对 相 界 面 变化 速率 的 影响 


Fig.4 Effect characteristics of supercooled effect on freezing 


rate of phase interface 


0.6} 


o4} 


Phase Interface (mm) 


L / 
27 
n 
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图 5 过 冷 度 对 凝固 相 界 面 变化 速率 的 影响 


Fig. 5 Effect characteristics of supercooling on freezing rate of 


phase interface 

图 6 显示 了 水 滴 凝固 过 程 不 同位 置 的 温度 变化 
| 线 。T01-T05 为 从 底面 到 表面 位 置 点 。 可 以 看 
出 ， 越 靠近 冷却 面 ， 温 度 降 低速 率 越 高 。 这 是 由 了 
凝固 第 二 阶段 主要 是 在 相 变 温度 下 由 热 扩 散 驱 动 而 
发 生 的 等 温 平衡 凝固 过 程 ， 导 热 是 该 过 程 的 主导 传 
热机 制 。 因 此 ， 越 靠近 冷却 面 ， 导 热 热 阻 相对 越 
小 ， 温 降 速率 也 相对 越 大 。 
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图 6 水 滴 不 同 高 度 位 置 的 温度 变化 特性 


Fig.6 Temperature variations of droplet freezing at different 


heights 

3 非 平衡 凝固 模型 
3.1 实验 系统 设计 

为 了 验证 所 发 展 模型 及 方法 的 有 效 性 ， 本 项 目 
自行 搭建 了 实验 系统 并 开展 了 相应 的 实验 研究 。 实 
验 系 统 示意 图 如 图 7 所 示 ， 由 半导体 制冷 系统 、 
Agilent 34970A 多 点 数据 采集 仪 、ANV TF100 PID 
温度 监视 器 、 温 度 控 制 器 、MotionXtra HG-100K 高 
速 摄像 机 、FLIR ”E60 红外 相机 、LED 无 影 光源 系 
统 、 工 控 机 及 调 压 电 源 组 成 。 实 验 时 由 半导体 制冷 
台 提 供 结 冰冷 环境 ， 由 PID 温度 控制 系统 将 制冷 台 
表面 温度 维持 在 实验 所 需 的 温度 条 件 ， 基 于 界面 追 
踪 的 方法 ， 采 用 高 速 摄像 机 记录 获得 水 滴 凝 固 过 程 
相 界 面 随时 间 的 变化 特性 ， 采 用 红外 相机 记录 水 滴 
凝固 过 程 温 度 的 实时 变化 特性 。 


Led Shadowless Lamp © 


mi 


High Speed Camera NOS y 
di N ^Supercooled Water 


Computer 


Temperature Monitor Temperature Controller 


图 7 水 滴 凝 固 特性 实验 系统 示意 图 
Fig.7 Experimental setup for droplet freezing 
32 ”实验 凝固 特性 与 计算 结果 的 比较 
在 水 滴 凝 固 实验 中 可 以 发 现 ， 在 过 冷 期 间 的 形 
核 与 枝 晶 生长 阶段 ， 水 滴 在 过 冷 条 件 的 驱动 下 逐步 
由 透明 态 向 模糊 态 过 渡 ， 当 水 滴 达 到 平衡 温度 ， 固 / 
液 界面 开始 出 现 并 由 液 相 区 向 固 相 区 移动 。 如 果 以 


fein B p HET 
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固 液 /界面 出 现 的 起 始 时 刻 作为 凝固 第 二 阶段 的 起 算 与 实验 相 界面 随时 间 的 移动 特性 吻合 较 好 ， 计 算 
点 ， 采 用 高 速 相机 记录 相 界 面 随时 间 的 变化 特征 ， 与 实验 温度 变化 历程 虽 存 在 一 定 偏差 但 变化 趋势 基 
并 同时 采用 红外 热 像 仪 记录 水 滴 不 同时 刻 的 温度 分 本 相似 ， 表 明 所 建 模型 的 有 效 性 。 

布 规律 。 图 8 显示 了 过 冷 度 为 -10C 条 件 下 凝固 第 二 12 
阶段 相 界面 随时 间 的 变化 过 程 ， 图 9 显示 了 相对 条 
件 下 相 界面 推进 过 程 水 滴 红 外 温度 场 随 时 间 的 变化 
特征 。 可 以 看 出 ， 在 凝固 第 二 阶段 ， 在 热 扩 散 的 作 
用 下 ， 相 界面 由 液 相 区 向 固 相 区 推进 ， 而 水 滴 的 低 
温 区 也 逐步 由 底部 向 顶部 发 展 。 
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图 10 实验 与 计算 凝固 速率 的 比较 


Fig. 10 Comparison of experimental and simulated freezing 


rates of phase interface 
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图 8 水 滴 凝 固 过 程 相 界 面 随时 间 的 变化 特征 
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Fig.8 Experimental evolution characteristics of phase interface 
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图 11 水 滴 凝 固 过 程 实验 与 计算 温度 变化 特性 的 比较 


Fig.11 Comparison of experimental and simulated temperature 


variations of droplet freezing 


4 结论 


本 文 针 对 飞机 结 冰 过 程 过 冷水 滴 结 冰 的 非 平衡 
凝固 效应 ， 开 展 了 过 冷水 滴 结 冰 凝 固 特性 的 计算 与 
实验 研究 。 主 要 结论 如 下 : 

图 9 水 滴 凝 固 过 程 温度 随时 间 的 变化 特征 (1) 基于 Enthalpy-porosity 模型 ， 发 展 了 过 冷 
Fig.9 Experimental temperature variations of droplet freezing 水 滴 结 冰 特 性 预测 模型 及 数值 计算 方法 ， 并 自行 搭 

图 10 显示 了 的 计算 与 实验 相 界 面 移动 特性 的 建 了 实验 台 ， 开 展 了 数值 方法 的 实验 验证 。 相 对 于 
比较 ， 图 11 则 显示 了 过 冷 度 分 别 为 -10'C 和 -14'C 条 传统 方法 ， 所 发 展 的 过 次 令 水 滴 结 冰 特 性 预测 模型 及 
件 下 ， 水 滴 表面 计算 温度 与 实验 红外 温度 变化 特性 mele EE 有效 表 征 非 平衡 凝固 效应 ， 因 而 可 用 于 


的 比较 。 从 图 10 可 以 看 出 ， 在 相同 过 冷 度 条 件 冷水 滴 凝 固 特性 的 预测 ; 
下 ， 改 进 预 测 模 型 的 计算 结果 与 实验 结果 的 吻合 程 (2) 在 结 冰 过 程 中 ， 过 冷 度 越 大 ， 水 滴 的 凝 
度 优 于 传统 T porosity 预测 模型 。 比 较 而 但 在 考虑 非 平衡 凝固 效应 的 条 件 


言 ， 考 虑 过 冷 效 应 的 改进 模型 预测 凝固 速率 要 高 于 下 ， 过 冷水 滴 凝 国 速率 要 高 于 不 考 谍 非 了 衡 凝固 效 
et UEDUD RU B e 吉 冰 过 程 中 过 冷水 滴 的 凝固 
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特性 具有 与 平衡 凝固 过 程 不 同 的 特征 。 有 效 表征 结 
冰 过 程 的 非 平衡 凝固 现象 是 改进 结 冰 热力 学 模型 ， 
提高 结 冰 特性 的 重要 基础 。 
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